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This work studied the actual and simulated technical performance between 
two grid-connected photovoltaic (GCPV) systems representing opposite 
latitudes. The system with a capacity of 5.4 kW> installed in Kelantan, 
Malaysia represents the northern equator, and the 183.6 kW» system installed 
in Cikarang, Indonesia, denotes the southern equator. The performance was 
simulated using PVsyst software, which included the energy output (Eoyt), 
reference yield (Y,), final yield (Yf), performance ratio (PR), and capacity 
factor (CF). The mean bias error (MBE) between the actual and simulated 
technical performance were as follows; for system A, the yearly MBE for the 
Eout: Yr, Yp, PR, and CF were -0.4%, 17.1%, -1.4%, -15.8%, and 1.4%, 
respectively, and for system B, the Egut, Y-, Y¢, PR, and CF values were 
9.80%, 18.3%, 10.0%, -7.2%, and 10.0% respectively. The results have 
proven that PVsyst has successfully simulated the yearly Egy, Y¢ and CF for 
both systems including PR, for system B, with MBE less than 10%. However, 
it is noteworthy to highlight that PVsyst significantly overestimated the Y, of 
both systems up to 18.3% and conversely underestimated the PR for system 
A by 15.8%, which highly likely caused by the Meteonorm imported weather 
data. 
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1. INTRODUCTION 


Renewable energy sources, such as solar, wind, hydro, geothermal, and biomass, are sustainable 
energy sources that have become popular in the modern world. Solar energy is one of the most essential and 
fundamental renewable energy sources because it is abundant and is used to generate electricity in most parts 
of the world [1]. Photovoltaic (PV) systems are practical and competitive options for transitioning to 
sustainable energy systems. Globally, the energy produced by PV systems has increased, reaching 760 GW in 
2020 [2]. The trend analyses show an exponential increase in the total installed capacity, along with increased 
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efficiency and decreased system prices, which piques the curiosity of more people globally and draws new 
stakeholders to a range of markets [3]. 

Researchers worldwide have studied the performance of PV systems from the module level to the 
utility-scale in various climates. The system requires International Electrotechnical Commission (IEC) 61724 
International Performance Monitoring to analyze the performance of several installations located and operating 
in various climates. Several technical performance parameters were recommended by the IEC through the 
IEC61724 standard to identify the best performers across the PV system. The performance parameters included 
the final yield ( Y¢), reference yield (Y,), and performance ratio (PR). These performance indices allow cross- 
comparison between several PV systems under different climatic conditions [4]. The technical performance 
evaluated using in-field monitoring data is also helpful for PV system designers, installers, researchers, and 
end users because they serve as a performance benchmark for component manufacturers [5]. The performance 
ratio for a 2.5 kW, installation in the Saharan environment in Algeria ranged from 66.66% to 85.93%. This 
broad interval is due to the significant temperature difference between seasons [6]. 

Anang et al. [7] performed a technical performance analysis of a 7.8 kW, grid-connected PV (GCPV) 
system rooftop in Kuala Terengganu, Malaysia. In addition, they conducted a performance analysis using the 
TEC61724 standard. According to the performance analysis, the best PR, value was 75.72%, and the average 
loss was 1.68 kWh/kW,/day. The annual capacity factor (CF) and overall system efficiency ranged from 13% 
to 16% and 10% to 12%, respectively. Furthermore, in Daher et al. [8], the PR for a 302.4 kW, installation 
operating in a dusty desert maritime climate ranged between 75% and 90%. The average daily array yield and 
final monthly yield were 5.1 kWh/kW, and 4.7 kWh/kW,, respectively. The average performance ratios for 
PV arrays and the global grid-connected system were 90% and 84%, respectively, corresponding to a monthly 
average of daily PV modules and system efficiencies of 12.68% and 11.75%. 

Lima et al. [9] reported another example of the standard application in the northeast region of Brazil, 
monitored from June 2013 to May 2014. The performance analysis of a 2.2 kW, PV system installed at the 
State University of Ceara, Fortaleza, revealed a suitable performance with an annual energy yield, an average 
daily reference, an array, and final yields of 1685.5 kWh/kW,, 5.6 kWh/kW,, 4.9 kWh/kW,, and 
4.6 kWh/kW, respectively. Moreover, a 600 W, PV system connected to the 220 V network of the Facultad de 
Ciencias Exactas building in Corrientes, Argentina, demonstrated a PR of 65% to 75% during the first ten 
months of operation between January 2011 and December 2012. This was consistent with those reported in 
other studies, specifically for systems installed on facades. However, the PR has fallen below 65% owing to 
network parameter instabilities since November 2011 [10]. 

The performance of a newly installed 281 kW, first GCPV solar farm in Lesotho, a country located 
in Southern Africa, was evaluated using IEC Standard 61724 parameters. The results demonstrated good 
performance, with a weighted PR of 70% compared with the global average of 70% to 80% for adequately 
performing PV farms [11]. It is critical to compare the performance parameters deduced from the actual data 
with the simulated performance to evaluate the accuracy of the software for future installations. Most 
frequently, the PV system performance using commercial software is performed during the project planning 
stage, particularly during the preparation stage [12]. Various PV simulation software packages employ various 
power models, databases, irradiance decomposition methods, and other features. Several simulation software 
packages, including HOMER, PV*SOL, RETScreen, and PVsyst, were used for designing PV systems. PVsyst 
is a design tool that optimizes grid-connected, standalone, and pumping PV systems based on their location on 
the map. The tool also depends on the consumer’s electricity profile and demand. It also provides financial 
visibility for the designed project, measuring environmental impacts in tons. The losses in the system can also 
be calculated [13]. 

Silva et al. [14] modelled a PV plant at the University of Campinas in Campinas, Brazil, using 
HOMER, PV*SOL, and PVsyst. According to their research, the value predicted by HOMER was under- 
predicted, PV*SOL was over-predicted, and PVsyst predicted the closest to the plant's measured values. These 
findings led to the recommendation of the PVsyst software for large-scale installations. Moreover, several 
studies have also applied PVsyst to estimate the GCPV system technical performance in equatorial regions. 
Tarigan et al. [15] simulated the techno-economic analysis of a GCPV system in Surabaya, Indonesia, using 
PVsyst and RETscreen. The technical analysis using PVsyst in the study shows that a 1 kW, GCPV system 
simulation managed to generate 1.3 MWh electricity per year, and the performance ratio for the system was 
approximately 73%. The study also found that the highest losses resulted from the array and inverter losses, 
which were 22% and 4.4%, respectively. 

Abdullah et al. [16] used the PVsyst software to perform a performance analysis of different types of 
PV modules for a 3 kW residential rooftop GCPV system in Selangor, Malaysia. Monocrystalline, 
polycrystalline and heterojunction with intrinsic thin layer (HIT) PV panels were used for comparative analysis 
in a 3.12 kW rooftop PV system. The PV module was analyzed using its exact location, weather, orientation, 
and losses to ensure a fair comparison. Polycrystalline produced the most energy to the grid, with 4046.9 kWh, 
followed by mono (3737.2 kWh) and HIT (3810.3 kWh). Although polycrystalline produced more solar energy 
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for the grid, the annual PR of the HIT panel was 81%, compared to 79% for polycrystalline and 79.5% for 
monocrystalline. The study concluded that the performance of HIT PV modules was superior to that of 
monocrystalline and polycrystalline -type PV modules owing to fewer losses and higher output [15]. 

In addition, Mansur ef al. [17] used PVsyst to design and simulate a 4.0 kW, solar PV system for a 
residential place under the Net Energy Metering (NEM) framework in Changlun, Malaysia. Based on the 
techno-economic study, the energy generated was 5704.4 kWh per year, with 9.7% used by the home load and 
the remaining 90.3% sent to the grid. The annual performance ratio was 79.6%, with an average daily energy 
production of 3.91 kWh/kW,. This NEM design setup is estimated to generate a profit of RM 1187 per year 
for residential customers while nearly reducing 4.0 tons of CO2 emissions into the atmosphere. 

Moreover, Hussain et al. [18] used the PVsyst software to perform a techno-economic analysis of 
commercial-size grid-connected rooftop solar PV systems under the NEM 3.0 scheme in Malaysia. At a net 
capacity factor (CF) of 18%, the system is expected to generate approximately 510 MWh of energy in the first 
year of operation. The annual degradation factors were 2.5% in the first year and 0.7% in the second year. The 
average PR was 0.803. The study found that temperature was responsible for the most significant losses, 
accounting for 7.77% of the total. In contrast, inverter losses were only 1.47%, which aligns with component 
manufacturing data This study also includes a financial analysis to assess the system's profitability, focusing 
on commercial buildings under the NEM 3.0 and resulting in an 8.4-year return on investment (ROI). 

The GCPV system generates the most significant amount of energy. Consequently, the performance 
of PV systems is a parameter that determines their efficiency. This study aimed to analyze the technical 
performance of two case studies of GCPV systems in the northern and southern hemispheres of the equatorial 
climate. The technical performance analysis was simulated using the PVsyst software and compared with the 
actual technical performance analysis as the benchmark. Nevertheless, this study will be limited to five 
technical performance parameters based on the data availability: Eoyt Y;, Y-, PR, and CF. Because most of the 
performance study was conducted within specific locations or geographical boundaries, this study was unique 
in comparing the technical performance of GCPV systems in different geographical boundaries yet different 
hemispheres of a tropical region, which are 6 degrees north and south of the equatorial line. 


2. METHOD 

The framework of this study is illustrated in Figure 1. The primary step was to compare the technical 
performances of the actual and simulated systems. This method framework will be explained in four sections, 
which are data compilation, PVsyst simulation, technical performance analysis using IEC61724:2021 and 
technical performances comparison for two GCPV sites. These two sites were selected because they are located 
at almost the same latitude but in different hemispheres. The systems were designated as systems A and B. 
System A is a5.39 kW, system in Kelantan, Malaysia, and system B is a 183.6 kW, system located in Cikarang, 
Indonesia. 

Section 2.2 discusses how these parameters were then used to simulate the system using the PVsyst 
software. To estimate the technical performance of the selected systems in the PVsyst software, the simulation 
design results must follow specifications according to the existing system: i) the model and specifications of 
the PV module and inverter used must be the same as the existing system; ii) the configuration of the array and 
inverter must be in line with the actual system; and iii) the configuration must fulfil the energy required for 
both systems. Then, in section 2.3, the technical performance of the existing systems is calculated based on 
TEC61724 standard, and lastly in section 2.4, the results of both methods are compared using error metric. All 
of these sections will be thoroughly discussed in the next section. 


2.1. Data compilation 

This section compiles the three types of information for systems A and B, as shown in Figure 2. The 
information compiled included system information, component information, and measured data from every 
5 min interval. The system information and component information will later be used in section 2.2, while the 
measured data will be used in section 2.3. 

System A is a 5.39 kW, GCPV system installed on the rooftop of a mosque in Kelantan, Malaysia. 
The system is located at 6.13°N and 102.26°E; it was installed by the end of 2014. The system consisted of 22 
silicon PV modules with 245 W, each. The 22 PV modules were arranged in two parallel strings, and the PV 
arrays were inclined at 10° facing the southwest. The strings are connected to PVI-5000/6000-TL-OUTD which 
is a 5-kW inverter which feed directly into the grid. 

Located at 6.29°S and 107.10°E, Cikarang, Indonesia, system B is a 183.6 kW, GCPV system on a 
factory's rooftop; it was installed by the end of 2020. The system consists of 540 silicon PV modules, each 
with a maximum power of 340 Wy. The 540 PV modules are arranged in 27 parallel strings, each consisting of 
20 PV modules. The PV arrays are inclined at 7°, facing south. Moreover, the strings are connected to Solid 
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Q-50 which is a 50-kW inverter which feed directly into the grid. A summary of the system information, PV 
module, and inverter specifications of both systems is provided in Tables 1-3, respectively. The measured data 
consisted of three parameters: in-plane irradiance (H;), module temperature (T,,), and AC power (P,,). The 
data were extracted from the data logger with 5 min logging intervals, in compliance with the requirements of 


the international standard IEC61724. 


Data Compilation 
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PVsyst Simulation 


Analysis using Photovoltaic system 
performance Part 1: Monitoring standard 
(IEC61724:2021) 


Comparison of actual and simulated 
technical performance. 
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Compilation of the GCPV system information of two case 
studies: 

A. 5,39 kWp system in Kelantan, Malaysia 

B. 183.6 kWp system in Cikarang, Indonesia 


Simulation of system A and B generated: 
i. Design result (Array and Inverter Configurations) 
ii. Technical Performance: (Eou,) (Y,), (Yd, (PR), and 
(CF) 


Actual technical performance: (Eou,) (Y,), (Yp, (PR), and 
(CF) were calculated. 


Mean bias error between the simulated and actual technical 
performance: (Eou,) (Yj, (Yp, (PR), and (CF) were 
calculated. 


Figure 1. Framework of the study 


System information: 
i. Coordinate 
Losses 
ili. Mounting type 
iv. PV module tilt angle 
PV module azimuth angle 
i. Array Configuration 
ii. Inverter Configurations 


. Commission Date 


Data Compilation: System A and System B 


Component information: 
i. PV module and Inverter i. In-plane Irradiance value 
model ii. AC power value 
PV module datasheet iii. Module Temperature 
iii. Inverter datasheet 


Measured Data: 


Figure 2. Data compilation of systems A and B (I continued to Figure 3) 


Table 1. General information for systems A and B 


Parameter/System A B 
Location Kelantan, Malaysia Cikarang, Indonesia 
Site coordinates 6.13°N, 102.26°E 6.29°S, 107.10°E 
Mounting Free - Standing Retrofitted 
PV array 5.39 183.6 
Array configuration 1x11 27X20 
Inverter configuration 1x1 1x3 
PV array tilt angle and orientation 10° facing southwest 7° facing south 
Commission date Dec-2014 Jan-2020 
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Table 2. PV module specification for systems A and B 
System A B 
Description Unit 


Brand Hanhwa Trina solar 
Model - Hanhwa- Q Cell Q-Pro-G3 TallMax TSM_DD4A (ID 
Type of PV technology Polycrystalline Polycrystalline 
Maximum power at STC Wp 245 340 

Open circuit voltage Vv 37.56 46.2 

Short circuit current A 8.85 9.50 
Maximum power voltage Vv 29.73 38.2 
Maximum power current A 8.32 8.90 
Temperature coefficient of voltage %°C!  - 0.33 - 0.29 
Temperature coefficient of current %°C! 0.04 0.05 
Temperature coefficient of power —-%°C! __-0.43 - 0.39 


Table 3. Inverter specification for systems A and B 


System A B 
Description Unit 
Brand ABB SMA 
Model - PVI-5000/6000-TL-OUTD Solid Q-50 
No. of inverter - 1 3 
Nominal power kw 5 50 
Maximum voltage Vv 600 1000 
Range of minimum and maximum MPPT voltage Vv 150..530 200..950 
Maximum AC A 25 80 
No. of MPPT - 2 3 
Inverter efficiency % 97 >98 


2.2. PVsyst simulation 

The PVsyst software was used in this study to predict the performance indices and energy output of 
systems A and B. The designs of both systems must be simulated before determining the performance indices. 
The steps applied for designing both systems using PVsyst are illustrated in Figure 3 [19]. The PVsyst 
simulation inputs are categorized into two sections: project information and simulation information. 

Project information consists of site details such as the project's name and coordinates. From the 
coordinates, the meteorological data of the selected systems were generated using the Meteonorm 8.0 database. 
The Meteonorm 8.0 database was selected based on a previous study that simulated the GCPV system using 
the PVsyst software [20]-[24]. 

Next, the simulation information consists of the systems’ main components, orientation, and losses. 
The PV module and inverter used are based on the installed system specifications in Tables 2 and 3. Because 
some of the PV module and inverter models are not listed in the PVsyst database, the specifications of these 
components were inserted manually into the PVsyst software. Then, the reference module temperature was 
applied such as most similar studies which a minimum of 20 °C and a maximum of 75 °C [25]-[27]. 

Then all these data were then simulated to obtain the design result. PVsyst processes the project 
information and the simulation information to generate the simulated GCPV design result, which comprises 
the PV array and inverter configurations. Then, the system losses, which consist of thermal losses, ohmic losses, 
light-induced degradation (LID) losses, soiling losses, incidence angle modifier (IAM) losses, auxiliary losses, 
aging losses, unavailability, and spectral collection, were inserted to enable the PVsyst software to simulate the 
design and performance results. The performance results included technical, economical, and environmental 
factors. However, owing to limitations in available information, this study is limited to analyzing the design result 
(array and inverter configuration) and technical performance, including Egyt.Y-¢, Y,, PR, and CF. 


2.2.1. Loss in PVsyst 

Various losses are involved in generating the energy output from a PV system. These include two 
types of losses: thermal capture and miscellaneous capture. Thermal capture losses are caused by cell 
temperatures higher than 25 °C. Miscellaneous capture losses are due to low irradiance, shading, dust 
accumulation on modules, mismatch, and wiring. In this study, the approach of determining the appropriate 
values of losses was prioritized using the values provided by the PV module and inverter manufacturer's 
datasheet, followed by values used in several works of literature in similar studies, and finally default values 
by the PVsyst software. The PVsyst simulation included 12 losses, as listed in Table 4. 

In this study, two losses were taken from the manufacturer's datasheet: the losses by aging (Kage) and 
the losses from the inverter (inv,,,,). System A has been operating since late 2014, approximately five years. 
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Thus, kage was simulated based on the power degradation information from the datasheet. However, system 
B has no losses due to aging because it is a newly installed system. For inv),);,, the values were obtained from 
the inverter datasheet based on the declared maximum efficiency. The other two losses obtained from the 
literature [28] are the losses due to dirt (Kgj-¢) and wiring (cablej9,,). Kgir~ was set to be 3% and cablejg.. 
was as also set to 3% based on [28] and many similar case studies conducted by [29]-[32]. 

For the other four losses, namely, temperature losses (Ktemp), mismatch losses (Kmm), LID, and IAM, 
the values applied were the default values provided by the PVsyst software. Notably, the thermal loss factor 
provided by PVsyst differs based on the mounting type of the array, which is 29 W/m°K for free-standing 
mounting and 20 W/m°K for retrofitted mounting. The temperature losses varied between the two systems, 
with system A experiencing an 8.86% loss as it was free-standing, while system B, which was retrofitted, had 
a 10.09% loss. The Kgpage loss was neglected in both cases because there was no shading on either selected 
site. Finally. The losses due to auxiliaries, unavailability, and spectral correction were not included because 


they were irrelevant in both cases. 


fo 


(4 
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Figure 1. PVsyst simulation steps 
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Table 4. Losses used in PVsyst and their references 


No PVsyst System A System B References 
Parameter Value Value Manufacturer datasheet PVsyst default Other 
1 Kage per year 0.60 0.55 
2 inVioss 3 3 V 
3 Kairt 3 3 V V [28] 
4 cablejocs 3 3 V [33] 
5° katie 8.86 10.09 V 
Thermal loss factor system A: 29 W/ m’K 
Thermal loss factor system B: 20 W/ m’K 
6 Kim 2.21 2.10 
7 LID 3 3 V 
8 IAM Losses 2.39 3.08 Vv 
9 


Kshade 0 0 V 
10 Auxiliaries - - : 
11 Unavailability - - 
12 Spectral correction - - 


2.3. Technical performance analysis 

The International Electrotechnical Commission (IEC) has developed parameters that have been used 
in various studies to analyze the performance of the GCPV system [34]—[38]. Several performance parameters 
are listed in IEC61724, but this study is limited to performing an analysis of (1) reference yield (Y,), (2) final 
yield (Y¢), (3) performance ratio (PR), and (4) capacity factor (CF) due to limitations in data availability. 

The reference yield (Y,) is the number of hours in which solar radiation must be at the reference 
irradiance levels to produce the same amount of incident solar energy observed during the reporting period 
when the utility grid or local load was available. The Y,. in (1) of a monofacial PV system can be calculated by 
dividing the total front-side in-plane irradiation by the module's reference plane-of-array irradiance, which is 
expressed as (1) [39]. 


— Hi 


Y.= 
r Gi 


(1) 

The H; is the total front-side in-plane irradiation, and G; is the reference plane-of-array irradiance, 
which equals 1 kW/m’. If the reporting period is equal to one day, then Y, would be, in effect, the equivalent 
number of sun hours at the reference irradiance per day. The Y; in (2) parameter was used to compare the energy 
output performance of PV systems in different geographical regions. This parameter is the ratio of the total 
energy generated by the system to the installed peak power of the PV array under rated conditions (STC), 
which is expressed as [39] (2). 


— Fout 
ears (2) 
Eout is the AC energy generated of the entire PV system per rated kilowatt system per year. Moreover, 
Eout is known as Eg,iq inside PVsyst. Next, Pp is the array power rating under the rated conditions (STC). 
Therefore, Yp normalizes the energy produced with respect to the system size. PR in (3), is the quotient of the 


system's final yield to its reference yield and indicates the overall effect of losses on the system, which is 
expressed as (3) [39]. 


Paes 

PR= “a (3) 

Y- is the final yield of the system and Y, is the reference yield of the system. PR represents how close 
to ideal performance a PV system is under real-world operating conditions; moreover, it allows for the 
comparison of PV systems, regardless of tilt angle, area, orientation, or nominal capacity. However, as the 
temperature of the PV module increased, most PR values were higher in colder months than in warmer months, 
resulting in additional losses. After accounting for energy losses, PR represents the amount of available energy 
[22]. PR typically ranges from 0.6 to 0.8, depending on the location, solar irradiation, and weather conditions. 

The annual CF in (4), is another performance parameter. This parameter represents the ratio of the 
annual energy output to the energy produced by the PV system if operated for 24 h per day at full rated power, 
which is expressed as (4). 


CF = —outveatly _ x 199 (4) 


Ppy ratedX365X24 
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CF is affected by the location of the PV system, and the higher the CF, the better is the performance of the PV 
system [37]. 


2.4. Result validation using mean bias error (MBE) 

MBE in (5) was the parameter used to assess the divergence between the simulated technical 
performance from the PVsyst and that measured from the GCPV system. Values near 0 are the best, negative 
values indicate the underestimation of simulated result from PVsyst and positive values indicate overestimation 
of simulated results [40]. It is expressed as (5). 

simulated value—measured value (5) 


MBE = x 100 


measured value 


3. RESULTS AND DISCUSSION 

This section is divided into three sections. Section 3.1 compares the actual and simulated irradiations 
for systems A and B. Then section 3.2 compares the actual and simulated technical performances of systems 
A and B. This comparison helps quantify the deviation of the simulated values from the actual values. Finally, 
section 3.3 is a comparison of the actual technical performance for systems A and B is presented to determine 
any significant difference in the technical performance for two GCPV systems located within the equatorial 
line but at different hemispheres. Nevertheless, it is very important to highlight that for system A, the available 
one-year actual data were for 2019. Meanwhile, the available one-year actual data for system B were for 2020. 


3.1. Comparison of the actual and simulated in-plane irradiation for systems A and B 

This study compares two sources of weather data: actual measured at the site and Meteonorm-derived 
long-term data (1985-2015) (PVsyst, 2022) for systems A and B. Figure 4(a) shows the comparison of monthly 
actual in-plane irradiation (Hj actuai) and simulated in-plane irradiation (Hj simulatea) Combined with the 
corresponding MBE for system A. The figure shows that the H; values of the actual and simulation for system A 
show a similar trend, although the magnitude of the MBE was quite high. The Meteonorm-derived data 
consistently overpredicted the value of actual irradiation throughout the year. The MBE varied in the range of 
11% to 32%. The highest MBE was in January, which was 31.9%. The highest H; recorded was in March for both 
Hj actual and Hi simulateas Which were 168.8 kWh/m? and 187.6 kWh/m* respectively. The lowest H; value was 
recorded in December for Hj actual 20d Hj simulateas Which were 103.4 kWh/m? and 118.0 kWh/m”, respectively. 

Figure 4(b) shows a comparison of Hj actual €0d Hj simulatea Combined with the corresponding MBE for 
system B. The trends of the actual and simulated H; were not in good agreement, except for three months: January, 
September, and November. The figure shows that Meteonorm-derived data also consistently overpredicted the 
value of actual irradiation throughout the year. The MBE varied in the range of 8% to 49%. The highest MBE 
was in February, which was 48.9%. The highest Hj actuay Was recorded in September with a value of 152.0 
kWh/m?, and in October for Hj simulated With a value of 168.8 kWh/m?. The lowest H; value was recorded in 
February for H; actuai and in January for Hj simulateas Which was 92.5 kWh/m? and 119.0 kWh/m”, respectively. 
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Figure 4. Comparison between the value of the measured in-plane irradiation (Hj gctuai) and the simulated 
in-plane irradiation (Hj simutatea) for (a) system A and (b) system B 


In summary, overprediction is expected from the Hj simutatea for both systems because Meteonorm 
provides long-term monthly data for a determined site created from the average statistical values over several 


Performance of grid-connected photovoltaic systems in Northern and ... (Yang Ilya Akila Abdul Rahim) 


866 0 ISSN: 2088-8694 


years. Thus, it will rarely be the same as the real values for any particular month, which also occurred in a previous 
study conducted by Vidal et al. in 2020 [23]. Although a significant value exists in the monthly MBE of the H;, 
the MBE for yearly H; for both systems is considerable as they range between 17% and 18%. The highest value 
Of Hi actuai recorded for system A is in line with the position of the sun since system A is located at the equatorial 
line, which is in March. Conversely, when the position of the sun is located further away from the equatorial line, 
which is in the tropic of Capricorn, Malaysians experience the northeast monsoon in November, which ends in 
February the following year; that is why the irradiation values in November and December are typically lower 
than those in the other months [41]. The irradiance ranged from 183 W/m? to 221 W/m2, which is lower than that 
of months such as March and April, which ranged from 244 W/m? to 259 W/m? [30]. 


3.2. Comparison of the actual and simulated technical performance results of systems A and B 

This section shows the comparison of the actual and simulated technical performances of systems A 
and B. The technical performance indices consist of Egy, Y-, Yp, PR, and CF which were presented on monthly 
and yearly basis. The monthly variations of the indices were illustrated in the form of a graph together with the 
corresponding MBE. 


3.2.1. AC energy generated 

Figure 5(a) shows a comparison of the monthly actual AC energy generated (Eout actual) and 
simulated AC energy generated (Eout simulated ) Combined with the corresponding MBE for system A. The 
figure shows that the actual and simulated E,,, trends are almost similar, with slight variations in January. The 
maximum Eo, was recorded in March for both Eout actual ANd Eout simulated » Which were 765.3 kWh and 737.2 
kWh, respectively, whereas the minimum E,y; was recorded in December for both Epout actuary and 
Eout.simulated » Which were 485.2 kWh and 474.6 kWh. The figure shows that the Eout simulated WaS very close 
to the actual values in the range of -0.4% to 5.6%, except for January, which was 11.9%. 

Figure 5(b) shows the comparison of the monthly Eput actual 22d Eout simulated COmbined with the 
corresponding MBE for system B. The actual and simulated E,,;¢ trends show that they are not in good 
agreement, except for three months, that is, January, September, and November. For system B, the maximum 
Eout actual WaS in September with a value of 22957 kWh, and the minimum was in February with a value of 
14550.63 kWh. It is observed that significant differences exist between the actual and simulated values, except 
for January, March, September, and November, with MBE less than 5%. The figure shows that, on average, 
the PVsyst software significantly overpredicted the Ej. simulated in the range of 0.2% to 29%, where the 
highest was in February. 
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Figure 5. Comparison between the value of the actual energy output AC energy (Eout actual) and the 
simulated energy output AC energy (Eout simulated ) for (a) system A and (b) system B 


The comparison shows that the maximum and minimum values of Eout actuai for both the systems are 
directly linked to the value of Hj actuai for each system. The highest Hj actual for system A was in March, and 
that for system B was in September, which is in line with the highest Eout actuai recorded. The same is true for 
the minimum value of Eout actual recorded. Therefore, the same conclusion can also be drawn for Eout simulated 
by PVsyst, as it also shows the same trend as the Eout actual. Although the MBE for monthly E,,; for both 
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systems is quite large, the yearly MBE only ranged between -0.4% to 9.8%. This MBE range can be considered 
small compared to a similar study conducted by Thotakura with a 30.64% deviation [42]. 


3.2.2. Reference yield 

Figure 6(a) shows the comparison of the monthly actual reference yield (Y; actual) and simulated 
reference yield (Y; simulated) Combined with the corresponding MBE for system A. The figure shows that the 
actual and simulated Y, values exhibit a similar trend. The highest Y, actuay Was recorded in March, with a value 
of 168.8 kWh/kW,, whereas the lowest was recorded in December, with a value of 103.4 kWh/kW),. The figure 
shows Y, simulated WaS Over-predicted by PVsyst, with the MBE varying in the range of 11% to 32%. The 
highest MBE was recorded in January with 32%. 

Figure 6(b) shows a monthly comparison of Y; actual 20d Y; simulatea Combined with the corresponding 
MBE for system B. The trends of the actual and simulated H; are not in good agreement, except for the three 
months: January, September, and November. The highest Y,. actual Was recorded in September (151.6 kWh/kW), 
while the lowest Y; actual Was recorded in February (92.5 kWh/kW,). The Figure shows that Y; simulated WaS 
overestimated for the entire year. The MBE varied in the range of 7.5% to 54%. The highest MBE was recorded 
in February, which was 54.3%. Analysing Y, for systems A and B, it is important to highlight that a high 
probability occurs that the simulated H; estimated using Meteonorm 8.0 in PVsyst is less accurate for both case 
studies conducted. 
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Figure 6. Comparison between the actual reference yield (Y; actual) with the simulated reference yield 
Y, simulated) for (a) system A and (b) system B 


3.2.3. Final yield 

Figure 7(a) shows the comparison of the monthly actual final yield (Y¢ actua) and simulated final yield 
(Y¢ simulatea)Combined with the corresponding MBE for system A. The figure shows that the trends of the 
actual and simulated E,,; are almost similar, with slight variations in January. The highest Y¢ actuay Was 
recorded in March with a value of 142.0 kWh/kW,, while the lowest Y¢ actuay Was recorded in December with 
a value of 90.1 kWh/kW,. The figure shows that Yr simulated predicted acceptable values for the entire year 
ranging from -0.3% to 10.9%. Nevertheless, note that the MBE in January was quite high, which was 10.9% 
that indicates the actual value is much lower than the simulated value. 

Figure 7(b) shows the comparison of monthly Y¢actua1 and Y¢ simulatea Combined with the 
corresponding MBE for system B. The actual and simulated E,,; trends show that they are not in good 
agreement, except for the three months, that is, January, September, and November. The highest Yr actua Was 
recorded in September, which was 119.5 kWh/kW,, while the lowest was recorded in February, with a value 
of 79.3 kWh/kW,. Fundamentally, Yp was directly proportional to the value of E,, produced in the same month 
[7]. Therefore, the higher the Hj actuai, the higher is the Eout actual and consequently, the higher is the Y¢ actual 
obtained. The figure also shows that the Y¢ simulated WAS Overestimated for the entire year. The MBE varied in 
the range of -0.2% to 34%. The highest MBE was in February, which was -33.9%. 


3.2.4. Performance ratio 

Figure 8(a) shows the comparison of the monthly actual performance ratio (PR actuai) and simulated 
performance ratio (PR simulated) Combined with the corresponding MBE for system A. The figure shows that 
the actual and simulated PR, trends were in good agreement with each other. The range of PR actual Was 
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84% to 89%, while the range of simulated PR, was 72% to 74%. The highest PR actual of 89% was obtained in 
October and the lowest PR actual Was in March at 84.1%. The MBE was in the range of -14% to -18%, which 
was in good agreement with a previous study conducted in the same climate region reported by 
Anang et al. [7], which was 15% to 16% 7. 

The highest MBE was attained in October, which was 18%. Because the highest Y¢ a-tua) Was recorded 
in March, it was expected that March will obtain the highest PR. Interestingly, the results showed that October 
was the month in which the GCPV system obtained the highest PR,, which indicates that Y¢ actual iS not the 
parameter that solely affects the PR value, as stated in a previous study [43]. The reason for the lowest PR value 
could be thermal losses, considering that March had the highest H; actual Value of 168.8 kWh/m? [23]. The 
average module temperature during the operation in the year 2019 ranged between 25°C and 67°C, and the 
average temperature of the PV module in March was 42 °C. Thus, it contributes significantly to the thermal 
losses experienced by system A. 

Figure 8(b) shows a monthly comparison of PR actua) and PR simutateaq Combined with the 
corresponding MBE for system B. The figure also shows that the actual and simulated PR, trends were in good 
agreement with each other. The range of PR actuai for system B lied between 77% and 86%, while the 
PR simulated ranged 74% to 75%. The highest PR actyay Was at 85.70% in February and the lowest PR actual Was 
in December at 77.24%. The figure also shows that PR gimulatea WaS lower than PR actuai- The MBE ranged 
from -3% to -13%. The highest MBE was in February, which was 13.1%. PVsyst underpredicted the value of 
PR simulated for the entire year in 2020. This may be because the of the overestimated value of the losses 
inserted into the PVsyst simulation. The value of the losses used was based on the worst-case scenario. 
Therefore, PR simulated 18 lower than PR actuai for system B. 
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Figure 7. Comparison between the actual final yield (Y¥¢ actua)) With the simulated final yield 
(Y¢ simulated) for (a) system A and (b) system B 
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Figure 8. Comparison between the actual performance ratio (PR actual) and simulated performance ratio 
(PR simulated) for (a) system A and (b) system B 
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3.2.5. Annual capacity factor (CF) 

The actual annual capacity factor CF,,tya; obtained in this study for system A was 15.7%, whereas the 
simulated capacity factor CF gimulatea Calculated from the simulation results with the Meteonorm-derived long- 
term average values was 15.5%. The MBE for the CF value was very small, 1.4%. The annual CF, ¢tya) obtained 
in this study for system B was 13.8%, whereas the CFimulateq WaS 15.4%, with an MBE of 10.1%. Because 
the value of CF is strongly dependent on the value of E,,,; recorded, as in (4), the results follow as the equation 
perceives. 

In summary, the ranges of the monthly actual and simulated technical performance are listed in 
Table 5. The table notably shows that for system A, PVsyst overpredicted the value of the H; and Y; in the 
range of approximately 11% to 32% but underpredicted the PR, by 15% to 18%. Similarly, for system B, 
PVsyst also overpredicted the value of the H; and Y; but in a bigger range of approximately 8% to 54 % and 
underpredicted the value of PR approximately 4% to 11%. 

Referring to Table 6, the yearly values of the simulated Eoyt, Y-, Yp, PR, and CF for system A were 
7.38 MWh, 1864.25 kWh/kW,, 1355.59 kWh/kW,, 72.8%, and 15.5%, respectively. In contrast, the actual 
yearly values of measured Ej, Y-, Y¢, PR, and CF were 7.40 MWh, 1591.6 kWh/kW,, 1374.7 kWh/kWp, 86.5% 
and 15.7%, respectively. The MBE for the Eout, Y-, Y, PR, and CF was -0.4%, 17.1%, -1.4%, -15.0%, and 1.4% 
respectively. 

For system B, the yearly values of simulated Eoyt, Y;, Ye, PR, and CF were 242.7 MWh, 1779.7 
kWh/kW,, 1325.43 kWh/kW,, 74.5%, and 15.1%, respectively. In contrast, the actual yearly values of the 
measured Eoyt, Yr, Ye, PR, and CF were 221.1 MWh, 1504.1 kWh/kW,, 1204.1 kWh/kW,, 80.3%, and 13.8%, 
respectively. The MBE for Eout, Y;, Y¢, PR, and CF was 9.80%, 18.3%, 10.0%, -7.2% and 10.0%, respectively. 

When comparing the monthly MBE with the yearly MBE of the technical performances for systems 
A and B, a larger deviation is observed for the monthly values compared to the yearly values. In the yearly 
overall analysis, unfavorable predictions in certain periods of the year may be masked by overestimations 
balancing out underestimations, as demonstrated in a prior study by Gonzalez-Pena [44] which found a smaller 
annual deviation ranging from 1% to 15%. 


Table 5. Range of monthly actual and simulated technical performance 
Range of monthly technical performance 

H; (kWh/m2) Eout (kWh) Y, (kWh/kW,)  Y; (kWh/kW,) PR (%) 
System A Actual 103.4 to 168.8 485.2 to 765.3 103.4 to 168.8 90.0 to 142.0 84.1 to 89.0 

Simulated 118.0 to 187.6 474.6 to 737.2 118.1 to 187.6 87.1 to 135.5 72.0 to 74.0 
MBE(%) Ss. to. 32.0 -0.4 to 11.9 11.1 to 32.0 -0.3 to 10.9 -15.0 to -18.0 
System B Actual 92.5 to 152.0  14550.6 to 21948.0 92.5to0 152.0 79.3 to 119.5 77.2 to 85.7 

Simulated 119.0 to 168.8 16329.0 to 22957.0 119.0to 169.0 89.0 to 125.0 74.1 to 74.9 
MBE (%) | 7.6 to 48.9 -0.2 to 29.4 7.5 to 54.3 -0.2 to 34.0 -3.9 to -10.5 


Table 6. Yearly actual and simulated technical performance 
i Yearly technical performance 
Eout (MWh) Y, (kWh/kW,)  Y; (kWh/KW,) PR (%) CF (%) 
7.40 86.5, 


System A Actual 1591.6 1374.7 15.7 
_ Simulated — 7.38 1864.3 1355.6 72.8 15.5, 

MBE (%) | i -0.4 17.1 -14 -15.8 -1.4 
System B Actual 221.1 1504.1 1204.1 80.3 13.8 
Simulated 242.7 1779.7 1325.4 74.5 15.1 

MBE (%) | 9.8 18.3 10.0 -7.2 10.0 


3.3. Comparison of the actual technical performance of system A and B with other studies 

In order to have the element of comparison, this section addresses the actual technical performance of 
systems A and B together with performances from other studies. The technical performance parameters of 
Y,, PR, and CF obtained in both systems; system A and B in this study, are found to be in the value range of the 
several studies as in Table 7. However, it is worthy to notify that both systems in this study have slightly lower 
CF when compared to the value of CF in tropical climates obtained by J. Ascencio- Vasquez, where the CF s 
ranged between 16% and 18% [45]. Despite the high irradiation levels in the tropical climate, the CF is lower 
than in other areas with high irradiation values owing to the year-round tropical rain, cloud cover, and constant 
high ambient temperature [45]. 
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Table 7. Technical performance of the previous studies including the present study 


Yearly technical performance 
No Author Location Climate Capacity (kW) Eout Y, Yr PR CF 


pus 
(MWh) (kWh/kW,) (kK Wh/kW,) (%) (%) 
1 17] Malaysia Af 78 9.4 - 1204.5 65.0 13.7 
78 10.7 - 1379.7 71.0 15.7 
2 [20] Malaysia Af 619.0 0.9 - - 77.0 12 
3 [46] Louisiana Cfa 1100.0 - - 1277.5 73.0 14.8 
4 [47] Malawi Cfa 830.0 - - 1551.2 79.5 17.7 
5 [48] Malaysia Af 18.0 9.9 - 1579.0 91.5 18 
18.0 8.9 - 1520.0 83.7 17.4 
18.0 8.9 - 1452.0 80.0 16.6 
6 [49] Ghana Aw 2500.0 3547.0 - - 70.4 16.2 
a [50] Thailand Af 3.5 - - 1387.0 76.4 - 
8 [51] Singapore Af 142.5 - - 1138.0 81.0 - 
9 [52] Malaysia Af 5.8 0.02 - 949.0 63.6 - 
10 Present study Malaysia Af 5.4 7A 1591.6 1374.7 86.5 1537 
11 Present Study Indonesia Af 183.6 221.1 1504.1 1204.1 80.3 13.8 


Legend: Cfa: Humid subtropical climate, Aw: Tropical Savanna Climate, Af: tropical rainforest climate 


4. CONCLUSION 

The study uniquely presented the performances from two GCPV systems located at opposites latitudes 
from the equatorial line; a 5.39 kW, free-standing GCPV system installed in Kelantan, Malaysia (northern 
latitude-system A), and a 183.6 kW, rooftop GCPV system installed in Cikarang, Indonesia (southern latitude- 
system B), which was monitored during 2019 and 2020, respectively. The study succeeded in comparing the 
simulated and actual technical performances of both systems. The following summaries and conclusions were 
drawn: i) For system A, the yearly values of simulated Egy, Y-, Ys, PR, and CF were 7.38 MWh,1864.25 
kWh/kW,, 1355.59 kWh/kW,,72.8%, and 15.5%, respectively; ii) In contrast, the actual yearly values of 
measured Eout, Y;, Y¢, PR, and CF were 7.40 MWh, 1591.6 kWh/kW,, 1374.7 kWh/kW,, 86.5%, and 15.7%, 
respectively; 111) The MBE for the Egy, Y;, Ys, PR, and CF was -0.4%, 17.1%, -1.4%, -15.8%, and 1.4%, 
respectively; iv) For system B, the yearly values of simulated Eout, Y;, Y¢, PR, and CF were 242.7 MWh, 1779.7 
kWh/kW,, 1325.43 kWh/kW,, 74.5%, and 15.1%, respectively; v) In contrast, the actual yearly values of 
measured Eoyt, Yr, Y¢, PR, and CF were 221.1 MWh, 1504.1 kWh/kW,, 1204.1 kWh/kW,, 80.3%, and 13.8%, 
respectively; vi) The MBE for the E,,;, Y,, Ys, PR, and CF was 9.80%, 18.3%, 10.0%, -7.2%, and 10.0%, 
respectively; vii) This study demonstrates that the most significant inaccuracy in the simulated technical 
performance was resulted from the value of irradiation imported into the PVsyst using Meteonorm 8.0; viii) 
Because Eoyt, Ys, and PR, are related to the amount of H;, some deviation occurs between the simulated and 
actual technical performance; and ix) While some difference exists between the technical performance of 
system A in the northern hemisphere and its simulated results, they are notably more consistent with each other 
than that in system B, in the southern hemisphere. 
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